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We calculate boundary operator product expansion coefficients C* r for boundary operators 0i b s c 
in the first column of the Kac table in conformal field theories, where {abc} specify the boundary 
conditions. For c = we give closed form expressions for all such coefficients. Then we general- 
ize to the augmented minimal models, giving explicit expressions for coefficients valid when 4>i t 2 
mediates a change from fixed to free (or any two specified) boundary conditions. These quantities 
are determined by computing an arbitrary four-point correlation function {4>i%<t>i% 1 <f>i d s 2 <j>i% 3 ) ■ Our 
calculation first determines the appropriate (non-logarithmic) conformal blocks by using standard 
null-vector methods. The behavior of these blocks under crossing symmetry then provides a general 
closed form expression for the C^„ r , as a product of ratios of F functions. This calculation was in- 
spired by the need for several of these coefficients in certain correlation function formulas for critical 
two-dimensional percolation and the augmented q = 2 and q = 3 state critical Potts models. 
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I. INTRODUCTION 

Conformal field theory (CFT) [1, 2] allows the determination of several results for 2D critical models in regions 
with boundaries. In the case of critical percolation these include crossing probabilities, such as Cardy's formula for 
the horizontal crossing probability [3, 4], Watt's formula for the horizontal- vertical crossing probability [5, 6], the 
expected number of horizontal crossing clusters [7, 8], and some new formulas for related crossing and other quantities 
[9, 10, 11]. Their derivations require the calculation of correlation functions of boundary condition changing (bec) 
operators 4>\ tS from the first column of the Kac table for the c = boundary CFT [3]. The resulting formulas include 
certain bec operator product expansion (OPE) coefficients C mnr , which (among other things) we calculate here. We 
first determine all C mnr of first-column Kac table operators 4>\ >s at c = (here there is no dependence on boundary 
conditions). As a byproduct, we find the (non-logarithmic) conformal blocks for an arbitrary four-point function 
(0i,2^ , i,s 1 ^ , i,s 2 0i,s 3 )- Then we generalize these results to the (augmented) minimal models, giving explicit expressions 
for the C^nr valid when <fii^ mediates a change from fixed to free boundary conditions. These are useful in formulas 
generalizing the exact factorizations for percolation found in [9, 11] to the Fortuin-Kasteleyn representation of the 
q = 2 and q = 3 state Potts models [12]. 

II. BOUNDARY OPERATOR PRODUCT EXPANSIONS COEFFICIENTS AT c = 

In this section we determine certain conformal blocks and operator product expansion coefficients at c = 0. These 
results are generalized to the augmented minimal models in section IV. 

In regions with a boundary there are three distinct classes of OPEs, namely bulk-bulk, boundary-boundary or bulk- 
boundary [13]. In this article we consider the boundary-boundary case exclusively. Correlation functions containing 
bec operators that belong to the Kac table obey differential equations, which allow us to determine them in closed 
form. From their solutions we can deduce the correlation functions and (our purpose here) the boundary OPE 
coefficients. As mentioned, our motivation is that these quantities arise in problems at the percolation point. 
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We begin with the tf>i t 2 Kac operator, which obeys the fusion rule 

[^1,2] X [(f>l tS \ = [<Pl,s-l] + [<Pl,8+l] ■ (1) 

Repeatedly applying (1) generates all Kac operators 4>i, s - This set of operators is closed under OPEs and constitutes 
the only elements of the Kac table that we consider. Note that (pi^ generates a "fixed to free" bcc in critical percolation 
[3], which is the origin of its interest for us. 

Since the (f>i s only depend on one parameter, we introduce the notation 

ipm '■= 01,m+l , (2) 

with m G N. The index m specifies the number of ipi fusions required to generate ip m . This includes the identity 
(f/'o)- In applications to percolation, the free-to-fixed operator is while 1^2 generates a cluster anchored at a point 
[9] and a physical interpretation of "03 is given in [10]. From the Kac weight formula we find that (at c = 0) the 
weights are given by 

m(m — 1) 1 f m 
h m = g 

Note that the integral values of h m are the Euler pentagonal numbers. 

The form of the OPE that we use in this paper (here including boundary condition dependence) is 

i+j 

iPf (A)Vf(O) = C$°\ h *- hi - h itp% c {0) +0(\) 1+h *- hi - h * . (4) 

p=\i-j\ 

p—\i—j\ even 

The sum includes all the primary operators dictated by the fusion rules for the ip m , and O(z) includes all contributions 
from the descendants of the primary fields. 

In a general theory the boundary OPEs depend on the boundary conditions assigned to each interval between 
the bcc operators. Thus there are additional indices on the OPE coefficients to represent the boundary conditions 
(see [13]). For percolation there are only two types of conformally invariant boundary conditions, free or fixed (see 
[3, 10]). However, they are equivalent by duality. As a result, the magnitudes of the OPE coefficients are independent 
of boundary conditions and their indices on the C mnr arc not needed. In addition, the normalizations of boundary 
operator two-point functions are independant of boundary conditions and we assume, with out loss of generality, that 
these two point functions are all normalized with coefficient unity, 

(^ m (x 2 ^ m (x 1 )} = (x 2 - xi)~ 2hm . 

(see sections III and IV for further details on this point). 
To extract the OPE coefficients, we make use of four-point correlation functions 

Km e hl (Vi(0)^(A)^(lM(0> = £ C.u.C^K H (A) . (5) 

V 

The functions fc; (A) are called conformal blocks. They are the contributions to the correlation function from 
individual conformal families appearing in the OPE (4), and are normalized so that the coefficient of the leading 
term (as A — > 0) is 1. The form of the correlation function in (5) follows on taking the limit A — > and using (4). 
Alternatively, one may take the limit A — > 1 and use the V'/c(l)V'j(A) OPE. The equivalence of these two expressions 
yields the crossing symmetry relation [1] 

p 9 

which we employ to determine the OPE coefficients. Note that the grouping of the indices in the conformal blocks is 
significant; however the indices of an OPE coefficient may be permuted without changing its value [1]. 

III. ANALYSIS 

We begin our analysis by calculating an arbitrary four-point function that includes at least one ipi , 

hm e hr (^m(0)M*)Ml)MO) = £<WW^l nrW ■ (7) 

P 
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Conformal covariance fixes the form of the four-point function up to a function of the cross-ratio [1]. We choose the 
form 



Gml,nr(x 1 ,X 2 ,X 3 ,X i ) = ^ m (x 1 )lp 1 (x 2 )lpn(x 3 )lpr(x i )) (8) 
= (a* - Xl )- h ™- h " +h r{x A - X^-^-^ixi - Xl )-hm+hn-hr F 



(x 3 - xi)(xi - x 2 ) , 
and pick x\ < x 2 < x 3 < £4, for later convenience. The tpi null- vector 

(3L 2 _ t - 2L_ 2 ) (9) 
implies that this function obeys the differential equation 

1 d \ / d N " ' 




" = ' 2 >' [(xi-xrf-^-xidxJ-'^dxJ I" 111,1 

where h\ = h m , h 3 = h n , and = h r . 

Exploiting conformal symmetry by letting {x\, x 2 , x 3l x^} — > {0, A, 1, 00} leads to 

° - F (A) + 3A(1-A) F (A) + 3A2(1 - Ap F(A) ■ (U) 

Note that our conditions on the Xi imply < A < 1. 

Equation (11) is linear, homogeneous, and has three regular singular points at 0, 1 and 00. Thus it can be 
transformed into a Riemann-Papperitz form with hypergeometric solutions [15]. We let 

F(X) = A m/3 (1 - A)" /3 /(A) (12) 

and insert the explicit forms of the weights (at c — 0), giving 

n _ , 2(l-2A + m(l-A)-nA) „ M , (m + n)(m + n + 1) + r(l - r) 

° " 7 (A) + 3A(T^ 7 (A) 9ATT^ /( } ' ( } 

(Note that up to (11) the analysis applies to the minimal models generally; only in (13) do we specialize to c = 0.) 
Defining 

m + n + r m + n + 1 — r 2m + 2 

a= , b= , and c= — - — , (14) 

we recover the hypergeometric differential equation 

= A(l - A)/" (A) + (c - (1 + a + b)X) /'(A) - abf(X) . (15) 
Using the two solutions of (15) and (12) gives the two conformal blocks 

m+n+r m+n+1 — r 



( m+n+r m+n+1 — r \ 
3 ajm+ii 3 Aj 



A and (16) 
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= A( 1 — )/ 3 (l - A)"/ 3 2 ^ ( 3 ^ 3 A J 



(17) 



When m = 2 (mod 3) the second hypergeometric solution (Eq. (17)) must be modified because of logarithmic 
effects, reflecting the logarithmic nature of CFT for c = 0. (Mathematically, this occurs because the parameter 
2 ( 2 ~ m } g — 2N.) However, this does not affect the validity of (16) since for blocks with weights differing by integers, 
the block with the highest-order leading term is regular [16]. (Four-point functions involving ip\ are also analyzed in 
[14], but this work does not consider the logarithmic case. ) 

Thus (16) and (17) give the (non-logarithmic) conformal blocks for an arbitrary four-point correlation function 
(with one ipi) °f the c = first column boundary operators ip s . (Explicit expressions for the null- vectors of these 
operators are given in [17].) 
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The crossing symmetry of conformal block (16) follows on using standard hypergeometric identities found, for 
example, in [18]. If we let m <-> n and A — * 1 — A in (16) and (17) the crossing symmetry is expressed in terms of the 
A — > 1 blocks as 

p ( 2m+2 \ p ( 2n-l \ p ( 2m+2 \ p / l-2n \ 

•r-m+1 M\ 1 V 3 / 1 V 3 / yn-1 (i _ \\ , 1 V 3 ) 1 \ 3 ) fn+l _ i\ /-i o\ 

• r ml,nr\ /K ) p ^ m + n + r ^ p ^ m _|_ n _|_i_ r \ - r nl, mrl 1 T p ^ rn+2-n-r ^ p ^ m+l-w+r ^ * r nl, mrl 1 'V ■ l ±c 7 

This allows comparison with (6), which in this case becomes 

Cl im ,m+lCr.n,m+lJ~™i^ lr {\) — C\^ n -.\^ n C r ^ l -\, m J : \ ll ^ mr (\ — A) — Ci i , i , n +iC r . n -|-i jrn ^ : '^ ? ^ I ,(l — A) . (19) 

Before proceeding further, note that for n = 2 (mod 3) there are logarithmic terms in the expansion of (16) around 
A = 1. In this case the crossing symmetry relation is 

p ( 2m+2 \ p ( 2n-l \ ?S W^ f m+l-n+r \ ( m+2-n-r \ 

rpm+l f\\ 1 I 3 / 1 V 3 / \m/3/i _ \\(l-n)/3 ^ 3 >q \ 3 >_q /-i _ \\g 

^tnl.nrl^ p / m+ti+r \ p f ro+n+l-r \ A V 1 'v | t A-2n \ V 1 'v 

V 3 / V 3 / q=0 V 3 /q 

+ 7 „ ^Vto .A m / 3 (1-A)"/ 3 V V 3 / q \ p- ^i(l-A)* (20) 

1 v 3 ^ 1 \ 3 y q =o q - \ q " r 3 r 

x(log(l~A)-V'(g+l)-V'(g+^)+V'(g+ g ^)+^(g+ m+ " 3 +1 ' r )) - 

with ^ the digamma function, instead of (18). 

By comparison with (4) we see that the exponent of the leading term in (20) belongs to the n — 1 block. Since 
the coefficient of this block is the same in (18) as in (20), we find that, regardless of logarithmic terms, the ratio of 
coefficients is 

n n ( 2m+2 \ p ( 2n-l \ 

v Ol,n-l,nW,ra-l,m L \ 3 ) L \ 3 ) / 91 \ 

ri-m,n,r ~ n r r / m+n+r\ r / m+n+l-r\ ' V ZJ J 

L'l.m^+l^r.n.m+l 1 ( 3 Jl ( 3 J 

There are four possible coefficient ratios that can be deduced from crossing symmetry, but (21) alone is sufficient to 
determine the OPE coefficients, so we omit the others. 
Now (21) gives 



1 ft p / m-\-n-\-r \ p / m+rH-l— r \ 

! ( i . n. in 4 I 1 I : M I >. J 



y/Tlm,n,r'R-n-l,m+l,r Cr,n-l,m /p ^ 2m+ l^ p ^ 2m+2 j p ( 2 "- X ) T ( 2 ^ ) 



a form which we will use to determine the OPE coefficients. 

We construct a general expression for C mnr , with < m < n < r, by using (22) to increase the indices in pairs. 
(Since we are free to permute the indices of C mnr , and the coefficients with an index are given by 

Comm = 1 , (23) 

this implies no loss of generality) It follows from (1) that m + n + r is even, and thus (m + n — r)/2 6 N. Now 

C^ )I2 +_ = 1. (24) 
Next, apply (22) to raise the second and third indices (m — n + f)/2 times. This results in 

m-n + T 



/p ^ ro+ra-r+2fc-l ^ p ^ 2k^ 
C m+n-r ^ m-n+r = ^ p ^ m + n-r + 2k j p ^ 2fc — 1 ^ ' (^^) 



Exchanging the first and second indices, and then applying (22) an additional (— m + n + r)/2 times, gives 



r in r(M=i)r(f) 



p p ^ m+n-r+2i ^ p / m-ra+r+2^-1 \ p ^ m+n-r+2e-l s 



/p ^ m+n-r+2fc-l ^ p ^2fe^ 
II \/ p ^ m+n-r+2k ^ p ^ 2fc^I J ) ' ( 26 ) 



k=l 
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Note that the order in which we used (22) to raise pairs of indices and obtain (26) is not unique; other orderings 
give different but equivalent expressions. 

Further manipulations of (26) lead to an expression that reflects the symmetry of the C mnr but is somewhat more 
complicated since it involves more products, 

r _ Ql V r gg)) (rgr r m) (rgr r (f )) (11 , r (f )) 



n-r 1 r(t))(n^ 1 r(D)(n^ 1 r i 



Equations (26) and (27) are the main results of this section. They provide closed expressions for all OPE coefficients 
between operators in the first column of the Kac table at c = 0. 

Considering particular indices can lead to more compact expressions. In the following, we find the lowest-order 
cases, i.e. C mnr with m — 1,2,3, directly from (22). 

Setting n = 1 and r = min (22) gives 



1 ' m ' m+1 " yr(M)r(i) ■ ^ 

Setting m — > m— 1 in (28) and permuting indices gives Ci im , m _i, which (as follows from (1)) exhausts the possibilities 
for Ci jm , r . If we let n = 2 and r = m + 1 in (22), include result (28), and redefine m — > m — 1 we find 



r(^i)r(l) r(l) 
2,m ' m ~ r(^) V r (t) r G)' {) 

If we let n = 2 and r = m — 1 in (22), redefine m — ► m + 1, and include result (28) we find 



^2, m ,n + 2 y (2m+1)r( M) ' 

Note that, by letting m —> m — 2, (30) also includes the coefficient C2.m-2.rn, which is the only other term in the 
fusion ip m il>2- The expressions in (31) and (32) below can be similarly manipulated to give all the non- vanishing C^mr- 
If we let n — 3 and r — m va (22) and include result (29) we find 



_ mr(f) / 2r(M) 

Gz.mjn+l — - W - ^ 2m+2 j~ ' (^^) 

If we let n = 3 and r = m — 2 in (22), redefine m — > m + 2, and include result (30) we find 



_ / 3(m+l) 

C3 ' m ' m+3 " V (2m+l)(2m + 3) ' (3 } 

Interestingly, this is a purely algebraic number. 

The results (28-32) also follow, as they must, from (26) or (27). The procedure leading to (28-32) can be generalized 
to obtain other coefficients in simpler forms than (26) or (27). 

The three coefficients used explicitly in [9, 10, 11] are 



r(|) _ 2 1 ' 2 ir 1 / 2 3 1 / 4 
;i)r(l " 

which follows immediately from (28) or (29) (and V function identities) 



Clu ~ ^ r m r'm ~ r(i/3)3/2 > (33) 



c r(|)r(|) _ 2^, 

Cl23 "V r(|) -3r(i/3)3/2- (34) 



from (28), and 



r _ r(f)r(§) / r(l) _ 2^ 
222 r(|) V r (l) r (l) 33/ 4 r(i/3)9/ 2 ' 1 j 



G 



from (29). 

The square of C112 appears as the leading coefficient in Cardy's formula [3]; our form agrees with Cardy's result. 
The ratio C112/C123 is a coefficient in a formula in [10]. The coefficient C222 appears as a proportionality constant 
in [9]. All of these coefficients are excellent agreement with simulations of the formulas in which they appear, for 
example C222 = 1-02993 ... as compared to 1.030 ± 0.001 from numerical results [9]. 



IV. GENERALIZATION TO MINIMAL MODELS 



In this section, by means of a straightforward generalization of section II, we calculate conformal blocks and operator 
product expansion coefficients for the (augmented) minimal models, i.e. we include coefficients for non-unitary minimal 
models. In the unitary case, many of the results that we derive are already known [13]. However, our extension to the 
augmented minimal models is to our knowledge new. As a consequence, some of our results are useful in calculating 
correlation functions of Fortuin-Kasteleyn clusters in q-state Potts models [12]. (For spin clusters in these cases the 
corresponding OPE coefficients may vanish.) In what follows, we present closed form expressions for the coefficients 
that are valid when (/> 12 mediates a change from fixed to free boundary conditions. (Actually, our results do not 
specifically require fixed or free boundary conditions. They also apply whenever only two boundary conditions, of 
any type, enter the problem and ipi^ specifies the points where they abut.) 

To avoid confusion with the index m used above, we denote the minimal model parameter by M, so that the central 
charge is 

C=1 -M(M+1) ■ (36) 
The boundary operators considered, as above, are generated by products of %p\ := 0i,2, namely 

nfnM - 2) 

Wn ■= <P\,n+U Wltn weights h n = M — . (37) 

Note that for the unitary minimal models the index n is restricted to < n < M. In this case M is an integer, but 
in fact our results apply for any real M with M > 2. 

We start with the four-point function, now including the boundary condition indices 

G a ^ d r a (xi,x 2 ,Xs, Xi ) := (CV)Vh b > 2 )C^3)Vv da (z4)> • (38) 

Proceeding as in section II, from (8) through (15), we find the conformal blocks (the boundary condition dependence 
of the correlation function enters only via the OPE coefficients, the blocks themselves being independent of boundary 
conditions) 



( M(l+m+n+r)-2 2(l+m+n-r) 
2(M+1) ' 2(M+1) 

M(m+1) 
M+l 

/ M(l-m+n+r) M(l — m+n-r)+2 

2 — m M n AT I — 

T m 7 X (X) = \2(M+T)(l - \) 2lifTT),Fi 2(M+1) ' 2(M+1) 

^ml.tirl^ A \ L ^) '2 r \ I 2+M(l-ro) 

M+l 



(39) 
(40) 



The Kac fusion rules impose the conditions 



m + n + r = 1 (mod 2) (41) 
n + r > m±l > \n — r\ . (42) 

For certain values of m and M, ^"^(A) is a logarithmic conformal block. This situation arises when (2 + M(l — 
m))/(M + 1) is a non-positive integer, so that (40) is undefined. However (39), which is the equation that we use, 
remains valid. 
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In the non-logarithmic case the crossing symmetry relations give 



( Mn-A 

-pm-L / \ \ _ V *" 1 ' / V M+1 / ^-ri-1 / 1 _ \ \ 



r ( 2-M(m-l) \ p 
! M+1 1 

-m+n+r) \ -p / M (l-m+ra-r)+2 
2(M+1) J 1 ^ 2(17+1) 

r / 2-Af(m-l) \ r / 1-Mn 

V M+1 J V M+1 / F »+i ri _,x 

, i+M(l-m-n-r) \ p f 2+M(l-m-n+r) ^ "LmrV ' 1 



2(M+1) / V 2(M+1) 

p / M(m+1) \ p / Mn-A 
-r-m+1 m = V M+1 J V Af+1 j t 



7m(i 



+m+w+r)-2 \ -p / M(\+m+n-r) 
2(M+1) / \ 2 (Af+1) 

M(m+1) \ p / l-Mn 



p / i^mti; \ p 

\ M+1 / V M+1 / 

' ' 2+M(l+m-n-r) \ p / A/(l + m-n+r) ^ nl ' mr ' ' ' 

1 ^ 2(M+1) 



2(M+1) 



As before we compare these crossing relations to the generic form 

/-iabc r~icda rv ac (\\ sidcb r~ibad rJ^^ *pn — l /l \\ ( AK\ 

b m,l,m+l u n,r,m+l a m+l- r ml,nr\ A / — ^n.l.n-l^m.r^-l^n-rti^mrl 1 A ) \*°) 

r~idcb /-ibad nfi^ 7T n +l /I 

^- y n,l,n+l^- , m,r,n+l^m+l*' nl,mrv / ' 

As explained above, except for percolation (M = 2) we must include the boundary condition labels on the OPE 
coefficients. We also need (as explained below) the a factors, which encode the normalization of the boundary two 
point function via 

(i;f(x 2 )^ b j a (x 1 )) = af(x 2 ~ Xl )- 2h > . (46) 

Thus the three point function is 

(^f(x 3 )^(x 2 )^( Xl )) = C^af{x 2 - Xl ) h >- h i- h *(x 3 - xtp-^-^ixa - x 2 ) hk ~ h ^ . (47) 
Symmetry under the exchange of operators in the three point function leads to the relations 

Cabc ac f~icba ac r^cab be /-ybac be r~tbca ab r^tacb ab {AQ\ 
ijk a k — L/ jik a k —^kij a j ~ ^ikj a j ~ u jki a i ~ Lj kji a i ■ V 10 ) 

When focusing on boundary operators, we could absorb the a coefficients into the normalization of the boundary 
operators as in section II. However the usual convention is to normalize the bulk operators to unity. Adopting that 
here keeps our results consistent with the literature [13], e.g. the unitary Ising model boundary OPE coefficients. 
Methods for evaluating the as are given in [13]. 
Thus we find 



r 



r 



/Af|m+!)\p (Mn=l\ dcb bad M 

\ M+i J L \ M+i J _ C n l n _ 1 C m r n _ 1 a n ^ 1 ^ 



' M(1 2y+5 r) ~ 2 ) r ( M( 2("M+lT r) ) C ^+X C ^,m + l<+l 



We focus on this particular relation because it is the only one of the four possible coefficient ratios that does not 
need to be modified in the logarithmic case. This can be shown by a series expansion analogous to that employed in 
section III. Thus (49) applies to all allowed m, n and r values. 

Using (48) and (49) we can write down the generalization of (22) for M ^ 2 



. . I r ( M(2+m+n+r)-2 \ r I M(2+m+n-r) 

#"£+l,n+l,r C^ l n+1 / agefc l { 2(M+T) ) l \ 2(M+1) 



JSabc f^abc \ n ad rv c< ^ 

m,n,r ^rn,n,r y m+1 n+1 



To simplify notation we have defined 



/p ( M(m+1) \ p / M(m+l)-l \ p / M (n+1) \ p / A/(n+l)-l \ 
V \ M+1 J L \ M+1 J L \ M+1 J L \ M+1 J 



(50) 



a 



jy-abc r~iabc j r 
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which is symmetric under simultaneous permutation of boundary conditions and operator labels, 



^bca 



ry-aoc _ jy-b 
m.n.r n.r.m 



K 



cab 



K 



cba 



K 



acb 



K 



hue 



(52) 

We can use (50) to raise pairs of operator indices on K^ c n r . 

From this point on, we consider only cases in which ipi mediates changes between free and fixed boundary conditions; 
thus all our boundary conditions are free or fixed. As a consequence, raising indices with (50) changes the type of 
intermediate boundary condition from one to the other. 

Using (50) on K^ c n r we can derive the relations 



Tfabc 
mnr 



K 



adc 



n 



v i M(r+£)-l \ p / Ml 
1 1 M+l I \ M+l 



\ 



-p / M (m-n+r+2l) \ p / 
1 \ 2(M+1) ) L \ 



M(m-n+r+2t)-2 
2(M+1) 



M{-m+n+r+2£) 
2(M+1) 



M(-m+n+r+2t)-2 
2(M+1) 



and 




v I M(-m+n+r+2k)-2 \ p / Mk 
1 1 2(M+1) J 1 \M+1 

p ,' M(-m+n+r+2k) \ p / Mfc-1 
' 2(M+1) J \ M + 1 



(53) 



(54) 



The fact that C§ ^ = 1 , for all ipj , a and b then implies that 



K, 



ddc 



c; 



ddc 



v dc 



and allows us to generalize (26) and (27). 
Thus we find 



mnr 



^m ^n 



V 



n 

i=i 



v i M(r+e)-l \ v ( Mi 
1 M+l ) \ M+l 



-p | M(m-n+r+2t) \ p / M{m-n+r+2l)-2 \ p / M (-m+n+r+2l)\ p / M (-m+n+r+2t)-2 
1 1 2(M+1) J 1 I 2(M+1) i 1 I 2(A/+1) J 1 I 2(M+1) 



n 

k=l 



M ( — m+n+r+2k) — 2 
2(M+1) 

'Ml — m+n+r+2fc) 
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A p / Mk \ \ 

\ p ( M(-m+n+r+2k) \ p / Mfc-1 A 
\ \ 2(M+1) y 1 ^ M+l j 

but (56) can once again be rearranged as 

« (nJT r (n^T r ( ^)) (ngT r (^)) (ng? 5 * r 



(55) 

\ 
/ 

(56) 



r^abc 
^ m.n.r 



J (rij=l^ ( M+l ) ^ ( M+l ) ) (1X7=1 r ^ -W- I" / ' ^.W-l. j J \y 1 ■-J 1 \ .W+l J ' \ .U+l 



(57) 



The question of whether the boundary condition d is free or fixed remains. By comparison with (55) we can see that 
if (—to + n + r)/2 is even (odd) than d represents the same (opposite) boundary state as c. For spin models, one can 
show that d is fixed (free) if there is (not) a Fortuin-Kasteleyn cluster that connects all three points. 
We list a few simple OPE coefficients that follow from (57), 



r<J J 



ff ff 

fFf V Q Q 2 _ FfF 



Ff 



°115 



a f f a FF 



a 



Ff 



\ 



r 



2M-1 
M+l 



r 



M 
M+l 



p / _2M_\ p / M-l 
1 ^M+l J 1 \^M+1 



C" 



fFF 
123 1 



x FF a fF 



- C Fff - 



12.! 



Jf„,f F 



\ 



( 3M-l \ p ( M \ 
\ M+l J 1 \^M+1 y 



3M 
M+l 



M-l 
M+l 



fff 
222 1 



FF 



a. 



./'./' 



FFF 
222 1 



FF 



\ 



p / 3M-1 
1 V M+l 



M 
M+l 



p / 2M A r Z' M-l \ p / 2M-1 
1 \^M+1 J 1 \^M+1 / \ M+l 



(58) 



(59) 



(60) 
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Here / (F) denotes free (fixed) boundary conditions. We do not need to differentiate between boundaries fixed 
in different spin states, since by symmetry the OPE coefficients must be independent of this. See [13] for more 
information on the as. One finds, for example, using the results of [13] for the necessary a values, for Ising clusters 
(M = 3), that C f 2 (i = 1.52552 ... and Cg£ F = 1.07871 .... 



V. CONCLUSIONS 



This paper derives closed form expressions for arbitrary Kac table boundary operator product expansion coefficients 
C mnr involving first-column operators <f>\ tS in the c = conformal field theory. Similar results for the coefficients 
Cmnr when 0i i2 mediates a change from fixed to free (or any two specified) boundary conditions in both unitary 
and augmented minimal models are also given. These coefficients occur in a variety of formulas for two-dimensional 
critical systems, derived via conformal field theory. These coefficients may also play a role in further exploring c = 
conformal field theory through descendant operator relations such as those in [10]. 
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